Propagation of an electromagnetic pulse through a waveguide with a barrier: A time 
domain solution within classical electrodynamics 
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An electromagnetic truncated Gaussian pulse propagates through a waveguide with piecewise 
different dielectric constants. The waveguide contains a barrier, namely a region of a lower dielectric 
constant compared to the neighboring regions. This set-up yields a purely imaginary wave vector 
in the region of the barrier ('electromagnetic tunneling'). We exactly calculate the time-dependent 
Green's function for a slightly simplified dispersion relation. In order to observe the plain tunneling 
effect we neglect the distortions caused by the wave guide in obtaining the transmitted pulse. The 
wave front of the pulse travels with the vacuum speed of light. Nevertheless, behind the barrier, the 
maximum of the transmitted pulse turns up at an earlier time than in the case without an barrier. 
This effect will be explained in terms of the energy flow across the barrier. The solutions obtained 
reproduce the shape of the pulses measured in the tunneling experiments of Enders and Nimtz 

PACS numbers: 03.50.De, 73.40.Gk, 03.40.Kf 



I. INTRODUCTION 



Tunneling is, like interference, a characteristic prop- 
erty of waves. This effect occurs both in non-relativistic 
quantum mechanics and in electrodynamics. Although 
the time dependent differential equations of both theories 
are fundamentally different in their structure, theoreti- 
cal calculations yield analogous results for the traversal 
time of the maximum of a wave packet's modulus. These 
calculations are based on the stationary phase approx- 
imation [jj^j as well as on a scattering ansatz ||. For 
sufficiently long barriers, the time delay is independent 
of the thickness and, thus, can correspond to arbitrary 
large effective velocities of the maximum of a pulse for 
crossing the barrier. This prediction is in agreement with 
results obtained by Nimtz and Enders in tunneling exper- 
iments with evanescent microwaves in a waveguide with 
a frequency below cut-off |^-^| . In these experiments the 
evanescent waveguide region, i.e. the barrier, is realized 
by an undersized region in between normal sized regions 
of a waveguide line. Due to the inhomogeneous cross 
section of the waveguide, we couldn't obtain an analytic 
expression for the transmission coefficient. Therefore, a 
recent microwave experiment || studied a barrier pro- 
duced by a low-dielectric-constant (£2) region which was 
placed in a rectangular waveguide of the same cross sec- 
tion, filled with a higher dielectric constant t\. This ex- 
perimental set-up is illustrated in Fig. [l]. The relation 
between wave number k and frequency lj is given by the 
dispersion formula (with the vacuum speed of light set to 



(1.1) 



0Jc/ e 2 possesses a real wave number k outside the bar- 
rier, but on entering the tunnel region with a lower di- 
electric constant k becomes imaginary, and the wave will 
spatially decay. 

In this paper we will consider the electrodynamic tun- 
neling for a barrier given by a variable dielectric constant. 
This set-up is amenable to a rigorous mathematical de- 
scription. It was theoretically investigated first by Martin 
and Landauer || . They concentrated on the tunneling of 
a Gaussian pulse with a narrow frequency range. Using 
a scattering ansatz they showed that the time delay of 
the center of mass depends only on the frequency deriva- 
tive of the phase of the transmission coefficient. For suf- 
ficiently long barriers, this delay becomes independent 
of thickness and thus corresponds to an arbitrary large 
effective velocity of the center of mass for crossing the 
tunnel region, which is known as the Hartman effect (lj] 
and has been experimentally verified first by Enders and 
Nimtz [El. 




where lo c represents the cut-off frequency of the empty 
waveguide and e is the variable dielectric constant pOj . 
Therefore a wave of frequency ujq with to^/ei < oj$< 



FIG. 1. The set-up of the considered waveguide. 



We want to examine how this effect is related to causal- 
ity. To carry out this goal, the fundamental solution of 
the given set-up, the retarded Green's function, will be 
constructed analytically assuming causality. This means 
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that this function vanishes outside the past light cone, 
i.e. the wave front travels with the vacuum speed of light. 
With the aid of this solution we will give an analytic ex- 
pression for the entire transmitted pulse, which has the 
observed superluminal property of the center of mass. In 
addition this allows us to determine the deformation of 
the pulse caused by crossing the barrier. It is assumed 
that the initial pulse is given by a truncated Gaussian 
wave packet located only to the left of the barrier. 



II. GREEN'S FUNCTION 



A. The Model 



The Green's function will be obtained using a La place 
transform. The structure of the dipersion formula (1.1) 
prevents an analytic inversion of this transformation be- 
cause of the different coefficients in front of oj 2 outside 
and inside the barrier, respectively. Therefore we con- 
sider in this paper a simplified model for the electromag- 
netic tunneling effect with dispersion formulas given by 



k = 



1' 



(2.1) 



with yfe coefficients dropped and cut-off frequencies 
mi = uj c /y/ei and TO2 = w c / outside and inside the 
barrier, respectively. The transmission coefficient of the 
barrier does not change qualitatively under this simpli- 
fication. The model given by these dispersion formulas, 
together with the Maxwell equations, still represents an 
electrodynamic case of tunneling. Within this model the 
velocity of a wave front limfc-»oo ui(k)/k = 1 is always 
given by the vacuum speed of light and it is assumed that 
the dielectric medium influences only the cut-off frequen- 
cies themselves. A similar model consisting of a classical 
scalar field which satisfies the one-dimensional relativis- 
tic Klein-Gordon equation with a rectangular potential 
barrier was investigated by Deutch and Low pd| ] . Under 
the condition that the tunneling amplitude is very small, 
they found an approximate solution given by a Gaussian 
wave packet, which turns up on the right hand side of 
the barrier as if its maximum took zero time to cross 
the barrier. Our goal here is to find an exact solution of 
the Maxwell equations yielding exact values for the tun- 
neling time of the maximum of a Gaussian wave packet, 
which can be compared to the approximate results given 
by Martin and Landauer. Furthermore, taking truncated 
wave packets of variable variance, the influence of these 
attributes on the tunneling time will be studied. 

Suppose that an electromagnetic pulse has been gen- 
erated in the region to the left of the barrier by an ap- 
propriate current which has vanished already before the 
wave front of the pulse reaches the left end of the barrier. 
The propagation of this pulse is then determined by the 
dispersion formulas together with the Maxwell equations 



with vanishing charge current. One can obtain the propa- 
gated field behind the barrier by solving these equations 
with the pulse located in front of the barrier as initial 
condition. This can be done with the aid of the retarded 
Green's function. In the case of the electromagnetic field, 
this fundamental solution can be written as an antisym- 
metric tensor Gxa/3, (a, = 0, . . . , 3), satisfying the four 
dimensional wave equation with the Dirac distribution at 
the fixed space-time position x as inhomogeneity 



(2.2) 



where □ = 8q — d\ — d 2 — <9f is the Laplace-Beltrami- 
operator of the Minkowskian space time p2| . Beyond 
that, G s must allow for the following boundary condi- 
tions: First, it has to vanish outside the past light cone 
at x, second its six components have to jump in the re- 
spective correct manner at the two boundaries between 
the propagating and the evanescent region and must take 
into account the metallic boundary conditions on the sur- 
face of the waveguide. In the next section, we will give a 
solution of this boundary value problem where the axial 
magnetic component G x 12 of Green's function shall serve 
as example. Given the relevant field components in front 
of the barrier at the time x° = 0, the expression 

B 3 {x) = J [G x i2,iE 2 ~ G xl2 , 2 Ei - G x i2, B 3 } dx 1 dx 2 dx 3 

N 

(2.3) 

yields the time evolution of the B 3 component of the mag- 
netic field behind the barrier x 3 > £/2 in terms of par- 
tial derivatives Gxap.-y of Green's function. The three- 
dimensional domain N = {x 6 R 4 |0 < x 1 < a, < 
x 2 < b,x° — 0} of integration is the interior of the whole 
waveguide line at time x° = 0. 



B. Analysis 

We solve the time dependent wave equation ( |2.2| ) by 
a Laplace transform relative to the time a; for a fixed 
x with the property a; > 0. The Laplace transform 
yields, after inversion, a function in the time domain 
which vanishes until a certain time is reached. Therefore 
this method is only suitable to get the advanced Green's 
function. The retarded Green's function, though, can 
then be obtained by a simple time reflection. 

The transformed components of Green's function are 
defined by 



G xa p{u) = £[G xa /3} — / e lLJX Gxapdx 



(2.4) 



where C denotes the Laplace transform with to in the 
upper half plane of complex numbers. 

For a waveguide with perfectly conducting walls given 
by the surface S, the boundary condition for the electro- 
magnetic field reads |l(J 
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(nxE) ls = 0, (n-B) 



\s 







(2.5) 



where n is the vector normal to the surface S. This con- 
dition is satisfied if one uses a Green's function which is 
invariant relative to the reflections at the four walls of 
the rectangular waveguide. This invariant Green's func- 
tion can be decomposed into the characteristic modes of 
the waveguide. The result for the Laplace transformed 
components can be written as 



i 00 

/it n h * ■ 

(2.6) 



t-,rj— — oo 



introducing the antisymmetric tensor 



/a/3 - 



ic\s 2 isic 2 ~sis 2 

-ic%S2 cic 2 ic\s 2 

-is\c 2 -cic 2 is\c 2 

s 1 s 2 —ic 1 s 2 —is\c 2 



with 



:= sin(x 1 t7r/a), 
:= sin(x 2 r]Tr / b) . 



(2.7) 



(2.8) 
(2.9) 



ci := cos(x 1 m/a), si 
c 2 := cos(x 2 rjjr I 'b) , s 2 

This tensor describes the x 1 - and x 2 -dependence of 
the modes represented by its respective pairs of num- 
bers (t, 77). Here we have introduced a one dimensional 
Green's function g a /3(x 3 ,x 3 ) to account for the depen- 
dence on the axial coordinate x 3 . Only this function 
contains the dependence on the frequency and will there- 
fore define the time evolutio n of the field. App lying the 
transformed wave equation (2J2) to Eq. ( |2.6| ), we find 
that the one dimensional Green's function has to be a 
solution of the ordinary differential equation given by 



d(x 3 ) 2 



+ lo 2 - m 2 (x 3 ) 



g a p{x 3 ,x 3 ) = 8{x 3 -x 3 ) 



(2.10) 

with the jumping cut-off freq uenc y determined by the 
simplified dispersion formulas (2.1), i.e. 



m{x 3 ) 



mi 
m 2 



u c /^, \x 3 \>t/2 
lu c /JT 2 , \x 3 \<t/2, 



(2.11) 



with uj c = ~K\J 1? /a 2 4- r/ 2 /b 2 for a rectangular waveguide. 
Notice that the cut-off frequencies mi and m 2 change 
with the type of mode. Thus the one dimensional G reen' s 
function also depends on the numbers (t, 77). In Eq. ( 2.10| ) 
no boundary terms at zero time emerge from the Laplace 
transform of the second time derivative due to the van- 
ishing of the advanced Green's function for a; < x and 
the property x° > 0. 



The one dimensional Green's function for (2.10) can 
be constructed from two linear independent solutions of 



the corresponding homogeneous equation. These can be 
chosen as a plane wave traveling from the left-hand side 
and the right-hand side, respectively, toward the barrier, 
i.e. 



Rape 
f- B af3 e 



■ik','J 



T, 



e 



ikx 




tape 



ikx' 



V ts a pe 

D „ikx 3 



x 3 < -1/2, 
\x 3 \<t/2, (2.12) 
x 3 > 1/2, 



x 3 < -1/2, 
\x 3 \<l/2, (2.13) 
x 3 > 1/2. 



with the total antisymmetric tensor e a p . The wave num- 
bers k and k are given by the dispersion formulas ( |2.l| ). 
Due to the upper half plane analyticity of a Laplace 
transformed function, the square roots in these equations 
have to be chosen analytic in the upper lu half plane. 
Hence the wave numbers always have a positive imagi- 
nary part. The respective coefficients R a p, A a p, B a p and 
T a p are uniquely determined by the boundary conditions 
at the two planes defining the evanescent region between 
x 3 = —t/2 and x 3 — 1/2. The coefficients of both solu- 
tions are assumed to be equal because of the symmetry 
of the barrier relative to x 3 = 0. Taking into account 
that the imaginary part of k is positive, the asymptotic 
behavior of the two solutions looks like 



lim 

c 3 — > — c 



,J a(3 



lim 

c 3 — *oc 



€l{x 3 ) = Q. (2.14) 



Therefore we can construct Green's function in the upper 
lu half plane as 



g a p(x 3 ,x 3 ) 



1 



W a t 



(2.15) 



with the Wronskian W a /3 of 



b {1) 

J af3 



and 



defined by 



W, 



a/3 



b i 2 Y_ 

a/3 dx 3 



a/3 dx 3 ■ 



(2.16) 



The case x 3 < x 3 with x 3 < -1/2 and x 3 > 1/2 is the 
interesting one for the investigation of the tunneling time 
problem for wave packets. In this case Green's function 
is suitable to calculate the field behind the barrier for a 
pulse assumed to be located initially only in the front of 
the barrier. In this case the Green's function reads 



g a/3 (x 3 ,x 3 ) = WTj-T a0 e 



ik(x 3 — x 3 ) 



(2.17) 



We now must consider the respective matching condi- 
tions for the components of the electromagnetic field to 
evaluate the coefficients of the two solutions. For simplic- 
ity let us only look at the calculations for the component 
£>3. This component has to be continuous at x 3 = —t/2 
and a; 3 = t/2. Furthermore we need the behavior of its 
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partial derivative perpendicular to the boundary-planes 
of the barrier. The transverse components B\ and Bi are 
continuous at any position of the two boundary-planes. 
Hence its partial derivatives parallel to the boundary- 
planes, i.e. the x 1 - and ^-derivatives, also have to be 
continuous at these planes. The equation divi? = shows 
that the x 3 -derivative of the component -B3 is continuous 
at the boundary-planes, too. These four continuity con- 
ditions together yield the coefficients of the correspond- 



(2). 



ing functions <fi 12 an d V12 



R 12 = D- X e- M (1 - e 2M ) (k 

B 12 = -D~ l 2k{k - K)e^ 3K - k W 2 

with the common denominator 

D = {k + K) 2 -{k~K) 2 e 2M . 



k 2 ) 



(2.18) 
(2.19) 
(2.20) 
(2.21) 

(2.22) 



Due to the absence of upper half plane zeros of D the 
coefficients are analytic in this half plane. 

We find the one dimensional retarded Green's function 
in the time domain after a time reflection. The inverse 
Laplace transform yields, for t :— x° — x° > 0, 



g a p{t,x ,x ) 



2tt 



g af 3(x 3 ,x 3 )duj (2.23) 



do -Ms 



with g a p(x , x ) given by Eq. (2.17 



In the case of the component B 3 we now h ave t o 
consider the transmission coefficient Tyi given in ( 2.21 ). 
Note that the expression k + k never vanishes in the up- 
per uj half plane. Thus, expanding the right-hand side of 
( [2.21D with l/(fc + k) 2 , we find 

T 12 = 77 i^_e^"^ . (2.24) 



(k + k) 2 



k — k 
k+K 



The term (k — k)/(Jc+ k) is bounded in the upper ui half 
plane by 1 and the imaginary part of k is always positive 
in this half plane. Thus the transmission coefficient can 
be written as a geometric series 

At this stage we need the simplified dispersion formulas 
( |2.l|) for an analytic inversion of the Laplace transform. 
Expanding each term of the series with (k + n) 2v cancels 
the ui 2 terms in the numerator. Thus the transformed 



Green's function (2.17) becomes 

~ gi 2(x 3 ,x 3 ) = -2iK 2 e^ s3 - x3 -^ 

1 (to 2 - m 2 ) 2l/ 
k {k + k) 4 "+ 2 ' 



V- 

z — t K - 



,ik£(2i/+1) 



(2.26) 



For our further calculations it is useful to write this as 
2 



012 (x ,x ) 



■£/■< 



l _ 



u=0 



(2.27) 



with the new function /„ given by 

/,M = (c 2 -m 2 ) e ^ 3 - 3 -' 



(2.28) 



Av+2 i(2v+\)i ^/uj 2 -ml 



iu+2 



and an effective cut-off frequency m = \J to 2 — to 2 . Dis- 
regarding the factor ahead of the fraction in ( 2.28[ ) for a 
moment, the function /„ corresponds to the time domain 
function [O 



j(x- -x-t 



= Q{t-{2v + l)£) 



t- (2u+l)£ 
t+(2v+l)£ 



2u+l 



■I 



iv+2 (moV< 2 -(2^+l) 2 ^ 2 ) , (2.29) 



with the step function Q(x) and the vth Bessel function J v (x). 

The additional exponential factor only causes a time translation. By performing the inverse Laplace transformation 
Or 1 one gets for q v (t, z) = Or 1 \] v (uS)(io 2 — m§) _1 ], with the abbreviation z — x 3 — x 3 > 0, the solution 

?„(*, z) = e(t-z- 2vl) ( ^ t ZlX\+% + +l)£ ) ^ Jiu+2 ( m oV(t- Z + i) 2 -(2iy+l) 2 £ 2 ) ■ (2.30) 

Because of the remaining factor lo 2 — tuq the second time derivative of q v (t,z) enters the function f v (t,z) = 
£ _1 [/„(w)]. Due to the vanishing of the function g„(0, z) for the values z + 2vl > no boundary terms arise from the 
Laplace transform of the second time derivative q^(t, z). Thus we have 



U(t, z) = -mlq v {t, z) - q v (t, z). 



(2-31) 



Now the inverse transform of the series in ( [2.27 ) yields for its terms, defined by h u (t, z) — C x \j v {\J up- — m 2 )], the 
expression [Of 
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h u (t, z) = f u (t, z) - mi I f v (V t 2 - u 2 , z^j Jx(miu)d% 



(2.32) 



Therefore we obtain for the integral ( 2.23| ) the series expansion 



*j 

9i2(t,z) = — VMM) 



v=0 



-mi 



q„(t,z) H o^(*' z ) 



q v {y/ t 2 — u 2 , z^j + -^kv {y/ 1 2 — u 2 , z^j 1 J\{miu)d 



(2.33) 



The structure of this Green's function can b e explained 
phenomenologically by looking at Eq. ( 2.32 ). Even the 
leading term q of the series, and thus ho, jumps from 
zero to a finite value at the boundary of the past light 
cone. This property of Green's function guarantees the 
causal propagation of every pulse. In opposition to the 
free space here the support of Green's function is not 
only the boundary of the light cone but its full interior. 
There are two reasons for this feature. First, consider the 



first term in Eq. (2.32) and accordingly the terms of the 
series q v given by Eq. ( 2.30 ). These terms contribute to 
the series only when t > z + 2vl. Thus the vth. term of 
the series represents the part of the field that leaves the 
barrier on the right hand side after 2v reflections at its 
boundaries. For the second reason take a look at the sec- 



ond term of Eq. (2.32). This term has nothing to do with 
the barrier itself but arises from the boundary conditions 
of the wave guide. Notice that the field is reflected there 
and back from the metallic boundaries while propagat- 
ing through the waveguide. This echo effect is described 
by the integral of the term in question. This integral 
represents a distortion in which all excitations that are 
noticeable at a given position in the time interval [0, t] 
take part. 



III. TUNNELING OF WAVE PACKETS 



A. Analysis 

Now we will solve the Maxwell equations for a given 
pulse using Green's function determined in the last sec- 
tion. Before starting with this calculation we want to 
simplify Green's function. It has been noti ced above that 
a distortion integral like that of Eq. ( 2.32 ) always arises 
in the case of guided waves. However, the tunneling ef- 
fect itself is described co mplet ely by the functions f v (t, z) 
in the first term of Eq. (2.32). We are interested in the 
undistorted delay induced by the tunneling effect only. 
Therefore, in the following calculations we consider only 
this term of Green's function, i.e. we set 



5i2 (M) 



00 r 1 



(t,z) 



(3.1) 



to determine the influence of the barrier itself on the 
transmitted pulse. With this simplification, the studied 
system has been mapped onto the model investigated by 
Deutch and Low . Their cut off frequency corresponds 
to our effective cut-off frequency mo = y m 2 — m\ . No- 
tice that, for comparing with experimental results, the 
contribution of the distortion term to the total transmit- 
ted pulse decreases with an increasing difference between 
the lower cut-off frequency mi and the central frequency 
Wo of the pulse. 

Using now Eq. ( |2.3| ), (2.6) with the simplified form 
T|) of Green's function with x° — 0, we obtain the 



component B3 behind the barrier for an arbitrary pulse 
started in front of the tunneling region. Our calculations 
will be carried out with a pulse given by a truncated 
Gaussian Hio-mode centered at some position x 3 = s < 
—£/2 in front of the barrier and with central frequency 
wq corresponding to the wavenumber ko 



In 



this case the relevant and non vanishing field components 
|10| are the real parts of 



E2(x) = iloq— sin(a; 1 7r/a) <p(x 3 ) 

7T 

Bz{x) = cos(x 1 7r/a) tp(x 3 ) 



(3.2) 
(3.3) 



with the Gaussian envelope 

^(x 3 ) = Q(- x 3 + s + 1 )e lkax \-~ {x3 ~ s)2/ ' j2 . (3.4) 

The distance between the maximum and the wavefront of 
the packet is given by the parameter 7 with the property 
< 7 < —s — 1/2. The upper boundary of 7 comes from 
the condition that the wavefront of the initial pulse has 
to be in front of the barrier at a; = 0. This is neces- 
sary because the barrier causes deformations of the pulse 
which are initially unknown. 

Wit h thi s initial pulse one can carry out the integration 
in Eq. (2.3) using the convolution theorem. Changing the 
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integration variable to v = t — x + £, we find the final barrier. With the relative coordinate u = x — x + s + £ 



solution for the component B3 of the pulse behind the one gets 

2 00 r 

s (x) = - cosl^w/a)^ cos(fc (w - s))T[}\u) + sm(k (u - s))T^(u) 
mi £ — ' L 



(3.5) 



with 



11+7 



(2i/+l)f 



2i/+l 



x J 4l /+ 2 (moV" 2 ~ 0+ !) 2 ^ 2 ) (cos(fc w)^(^ 2 )(M - u) T sin(/c w)^( 2 ' 1 ) (u - rfw, 



and the abbreviations 



(3.6) 



2 / 4i 2 6 

i?W (t) = t wg - m\ + 2fc (fc + wb) - —r + -7 



^ (2) (t) = fco (c 2 - 



2 12t 2 6 \ / 2 2 4£ 2 2 

m 2 - — + ^ J I ^0 - ™ 2 - ^4 + ^2 / 



(3.7) 
(3.8) 



The integral in Eq. ( |3.6| ) can be evaluated only numeri- 
cally. Note that, for a fixed value of u, all of the functions 
(1 21 

Vli ' (u) with an index v larger than some index vq van- 
ish. Naturally, the index vq depends on the value of u. 
So we have to do only a finite number of numerical in- 
tegrations to obtain the component B3. The envelope of 
B3 is shown in Fig. || for different barrier thicknesses £ 
and fixed central and cut-off frequencies. The distance 
7 between the maximum and the wave front of the wave 
packet has to be large enough to prevent deformations 
of the transmitted pulse that arise in the case of a pulse 
with too large high-frequency components. For this and 
the following results 7 has been chosen to be five times 
the initial variance <7q of the pulse. Due to the trivial 
dependence of the solution on the coordinate x 1 we have 
always set x 1 — corresponding to the boundary of the 
waveguide. Because we have eliminated the echo effect 
caused by the waveguide itself, the pulse does not change 
its shape outside the barrier region. Thus the graphs in 
Fig. correspond to the time evolution of the wave packet 
measured by an observer at an arbitrary position behind 
the barrier. To determine the tunneling time of the max- 
imum of the wave packet, we can use the free propagation 
of the wave packet outside the barrier, i.e. the fact that 
it travels there with the vacuum speed of light without 
changing its shape. Then the maximum of the packet ar- 
rives at the left end of the barrier at a; = —s — £/2. Now 
let r be the value of u at which the envelope of B^(u) 
has its maximum. Considering the definition of the co- 
ordinate u, one obtains for the arrival time at the right 
end of the barrier x° = t — s — £/2. Thus the tunneling 
time is given by r. 
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FIG. 2. Graphs of the envelope of the component B3 of the 
magnetic field for Gaussian wave packets transmitted below 
cut-off across barriers of different thickness £ as a function of 
the relative coordinate u. The cut-off frequencies are mi = 1 
and m,2 = 4. In all cases coo = 3.2 and a — 10. The graphs 
are scaled by the factors in the brackets. 



The graphs in Fig. g show that the transmitted pulses 
are also Gaussian-like wave packets, but exponentially 
damped with growing barrier thickness. To obtain the 
tunneling time and the variance of the packets we have 
fitted Gaussian wave packets to the graphs of Fig. ^. The 
resulting values for r and the variance a are listed in 
Table I. Furthermore, this table shows values also corre- 
sponding to Gaussian-like solutions for other cut-off frc- 



G 



quencies of which the graphs are not shown here. While 
in the case of i = 0.5 the tunneling time corresponds to a 
subluminal average velocity of the maximum, it increases 
more and more slowly with growing barrier thickness. 
This result agrees with the experimental observation that 
for sufficiently long barriers the tunneling time is inde- 
pendent of thickness || ■ The variance of the transmitted 
wave packets decreases with increasing barrier thickness. 

To determine the dependence of the tunneling time 
on the central frequency of the initial pulses we have 
calculated the corresponding transmitted wave packets 
for a barrier of fixed thickness I — 5 with m\ = 1 and 
ni2 = 4. We have considered Gaussians with variances 
both above (<7o = 10) and below (<jq = 4) the thick- 
ness. The resulting pulses are also Gaussian-like except 
for the low- variance-pulses with ojo ~ 7712/2. The param- 
eters of the Gaussian solutions are given by Table II. The 
tunneling time increases with the central frequency but 
corresponds always to a superluminal average velocity for 
the maximum of the pulse. The maxima are shifted to 
higher values of r with increasing loq because of contribu- 
tions from Fourier components above the barrier cut-off. 
These components are also responsible for the slightly 
higher times of the narrow packets and its distortion at 
higher central frequencies. The variance of the transmit- 
ted packets decreases with increasing uj . 

TABLE I. Numerical results for tunneling times and vari- 
ances of wave packets that have crossed barriers of different 
length. Two different pairs of cut-offs are considered with a 
corresponding central frequency in between. In both cases the 
initial variance is <xo = 10. 



account for this fact, we compare the tunneling times of 
our echo-free solutions with that of Martin and Landauer 
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Let us now compare the tunneling times obtained from 
our solutions with that of Martin and Landauer ||. 
These authors have pointed out that the time delay of 
the center of mass for a pulse restricted to a wide vari- 
ance in the time domain depends only on the frequency 
derivative of the phase a of the transmission coefficient, 
i.e. 



1"ML 



da 



,dk 



(3.9) 



where the expression has to be evaluated at uj = luq. Due 
to the symmetry of the transmitted pulse the time delay 
for the center of mass and the maximum are the same. 
But the time tml does not represent the pure time delay 
of the maximum caused by the barrier itself: Due to the 
echo effect of the waveguide which also affects the trans- 
mission coefficient, tml is shifted to higher values. To 



by taking in Eq. (3.9) the phase of the transmission coef- 
ficient with mi = 0, corresponding to a free propagation 
outside the barrier. The graph of tml and our values of 
r are plotted in Figs. || and |j as functions of the central 
frequency and the barrier thickness, respectively, for the 
parameters considered above. A difference between the 
two tunneling times arises only for higher central frequen- 
cies because of the growing contribution of high Fourier 
components to the transmitted pulse. Notice that the ap- 
proximate result of Martin and Landauer is valid only for 
pulses with a narrow frequency range. Furthermore our 
values for the tunneling time retains a small dependence 
on thickness, but corresponds nevertheless to a superlu- 
minal velocity for the maximum of the pulse. At this 
moment it should be emphasized again that the wave 
front travels always with the vacuum speed of light. Be- 
low, the possibility of superluminal maxima within an 
underlying causal propagation will be explained in terms 
of the energy flow across the barrier. 

TABLE II. Numerical results for tunneling times and vari- 
ances of wave packets with different central frequencies be- 
tween both cut-offs and initial variance <jq = 10 and <7o = 4, 
respectively. The barrier thickness is I = 5. 
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FIG. 3. Graph of the tunneling time for a wave packet 
calculated by Martin and Landauer || from the frequency 
derivative of the phase of the transmission coefficient as a 
function of the central frequency of the packet. Solid line: 
The cut-off frequencies are mi = 1, mi = 4 and the barrier 
thickness is I = 5. The dashed line corresponds to a vanish- 
ing cut-off frequency outside the barrier (mi = 0). The dots 
represent the tun neling time of the maximum of our solution 
given by Eq. (3.5) with ctq = 10. 
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FIG. 4. Graphs of the tunneling time as given both by |j] 
with vanishing cut-off frequency mi out side the barrier (lines) 
and by the maximum of the solution (3.5) (dots) as a func- 
tion of the barrier thickness for two different pairs of cut-off 
frequencies and a central frequency in between. The dots 
correspond to the values given in Table I. 



B. Interpretation 

Apparently, the maxima of the solutions we have given 
above cross the barrier with a superluminal velocity. In 
this section, we want to explain this property of the trans- 
mitted wave packet in terms of the energy flow across 
the barrier. In Minkowskian space-time it is given by the 
integral curves of the 4-vector field T Q with the electro- 
magnetic energy density T° = j{E ■ D + B ■ H) and the 
Poynting vector T j = (E x H)j for j = 1, 2, 3. 

The fact that T Q is a continuous differentiable vector 
field induces an important property of its integral curves 
in space-time: Curves with different initial positions do 
not intersect each other. This property is not destroyed 
by the jumping conditions for the electromagnetic field 
at the ends of the barrier. To prove this claim we have 
to show only that any curve flowing on a boundary of 
the barrier from one side has a unique continuation at 
the other side of the boundary. This means that there 
are no branching points for the curves at the boundary. 
Due to the continuity of the component T 3 at the bound- 
ary these points could only arise if one assumed that the 
component T 3 vanished there and thus on both sides 
the curves were tangential to the boundary at this posi- 
tion of space-time. This would be the case for a curve 
that comes, for example, from the right, is tangent to 
the boundary and than goes back to the right. At the 
position where the curve is tangential to the boundary, 
another curve could flow into this curve from the other 
side of the boundary, leading to a branching point. But 
this situation is not possible because the tangential vec- 
tors of the curves have to point in the directions of the 
future light cone. Looking at the neighboring curves of 
those considered above, these possible directions would 



be inconsistent with the continuity of T 3 and the fact 
that the curves do not intersect each other outside the 
boundary. Thus there are no branching points at the 
boundary. 



bifurcation curve 




wave packet 



barrier 



FIG. 5. A qualitative plot of the energy flow in a 
x°-x 3 -plane of the Minkowskian space-time. The plot in- 
tensity of the initial wave packet corresponds to its energy 
density. 



So we obtain the qualitative picture for the integral 
curves of the energy flow in the space-time shown in Fig. 
[| The curves originate in the initial pulse of which the 
wave front is still to the left of the barrier at time x° = 0. 
The curve originating in the wave front of the pulse has a 
slope of one because the wave front propagates with the 
vacuum speed of light. The fact that the curves do not 
intersect each other allows the initial pulse to be decom- 
posed into two connected parts from which transmitted 
and reflected curves, respectively, originate. If there is 
no energy absorbed by the barrier, the starting position 
s of the bifurcation curve separating transmitted and re- 
flected curves is given implicitly by the transmitted por- 
tion of the energy of the initial pulse, i.e. 



1 \x° 



(3.10) 



The integration along the time axis can be done at any 
position x 3 = p behind the barrier. 

Now, by means of the energy flow as shown in Fig. [|, 
we want to explain the existence of solutions with su- 
perluminal maxima within a causal theory. The part of 
the initial pulse between its wave front and the starting 
point s of the bifurcation curve is mapped along the en- 
ergy flow on the time axis at x 3 = p. The closer two 
neighboring points of the initial pulse at x° — are to 
the starting point s of the bifurcation curve, the more 
the distance between them is extended by this mapping. 
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This is necessary because of the arrival of transmitted 
curves at x 3 — p also at arbitrary late times. Due to 
this spreading of the curves the energy density of the 
transmitted pulse at x 3 — p begins to decrease from a 
particular time corresponding to the arrival time of the 
maximum behind the barrier. In other words, the trans- 
mitted pulse results from a redistribution of the energy 
contained in the forward tail between the front and s of 
the starting pulse. Thus the maxima of the initial and 
transmitted pulse are not causally related. 

The whole picture arises as a consequence of the math- 
ematical claim that the curves do not intersect each 
other. These curves themselves, of course, cannot be ob- 
served in any experiment. But we believe that they are a 
suitable tool to get a classical picture of the mechanism 
of the tunneling effect. Within this classical interpreta- 
tion the surprising result of our solutions is the almost 
exact reconstruction of a Gaussian wave packet behind 
the barrier. This effect yields a pulse-reshaping 0]. 

To obtain a more physical point of view one can ask 
at which time the transmitted pulse exceeds an arbitrary 
threshold behind the barrier. Due to the damping and 
squeezing of the transmitted pulse this happens always 
at a later time compared to a pulse which crosses no bar- 
rier. That means for an observer behind the barrier that 
he would not detect the tunneled pulse earlier than the 
freely propagated one, in agreement with causality. 



Within this interpretation the energy of the transmit- 
ted pulse can only originate from a connected part be- 
hind the wave front of the initial wave packet because the 
integral curves of the energy flow do not intersect each 
other. In this sense, the Gaussian shape of the transmit- 
ted pulse can only be interpreted as an amazing interfer- 
ence effect. Due to the propagation of the wave front with 
the vacuum speed of light, it is not possible to obtain a 
superluminal maximum if the barrier thickness exceeds 
the distance between the maximum and the wave front. 
Thus, in the case of truncated wave packets, the time 
delay of the maximum does not stay independent of the 
thickness for all barrier length. In summary, the results 
of the tunneling experiments can be obtained from the 
causal Maxwell theory. 
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IV. SUMMARY AND CONCLUSION 

In Eq. ([5^), we have given an exact analytic expres- 
sion in the time domain for the causal Green's function 
of a model that describes an ideal case of electromag- 
netic tunneling. The structure of this function allows for 
a reduction to those terms which describe only the pure 
tunneling effect without the distortions caused by the 
waveguide. With this reduced Green's function, we calcu- 
lated the shape of transmitted wave packets for truncated 
Gaussians as initial pulses. The resulting pulse can be 
also Gaussian-like. In agreement with the approximate 
result of Martin and Landauer Q and the experiments 
on microwaves by Enders and Nimtz [^| , the time delay of 
the maximum of the pulse becomes nearly independent of 
the thickness for sufficiently long barriers. Furthermore, 
the variance of the transmitted packet decreases with in- 
creasing barrier thickness. By examining the properties 
of the energy flow, we have found consistency of a super- 
luminal pulse maximum with the causality of Maxwell's 
theory. 
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